home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / slaed6.z / slaed6
Encoding:
Text File  |  1998-10-30  |  2.7 KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSLLLLAAAAEEEEDDDD6666((((3333FFFF))))                                                          SSSSLLLLAAAAEEEEDDDD6666((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SLAED6 - compute the positive or negative root (closest to the origin) of
  10.      z(1) z(2) z(3) f(x) = rho + --------- + ---------- + ---------  d(1)-x
  11.      d(2)-x d(3)-x  It is assumed that   if ORGATI = .true
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE SLAED6( KNITER, ORGATI, RHO, D, Z, FINIT, TAU, INFO )
  15.  
  16.          LOGICAL        ORGATI
  17.  
  18.          INTEGER        INFO, KNITER
  19.  
  20.          REAL           FINIT, RHO, TAU
  21.  
  22.          REAL           D( 3 ), Z( 3 )
  23.  
  24. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  25.      SLAED6 computes the positive or negative root (closest to the origin) of
  26.                       z(1)        z(2)        z(3) f(x) =   rho + --------- +
  27.      ---------- + ---------
  28.                      d(1)-x      d(2)-x      d(3)-x
  29.            otherwise it is between d(1) and d(2)
  30.  
  31.      This routine will be called by SLAED4 when necessary. In most cases, the
  32.      root sought is the smallest in magnitude, though it might not be in some
  33.      extremely rare situations.
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      KNITER       (input) INTEGER
  38.                   Refer to SLAED4 for its significance.
  39.  
  40.      ORGATI       (input) LOGICAL
  41.                   If ORGATI is true, the needed root is between d(2) and d(3);
  42.                   otherwise it is between d(1) and d(2).  See SLAED4 for
  43.                   further details.
  44.  
  45.      RHO          (input) REAL
  46.                   Refer to the equation f(x) above.
  47.  
  48.      D            (input) REAL array, dimension (3)
  49.                   D satisfies d(1) < d(2) < d(3).
  50.  
  51.      Z            (input) REAL array, dimension (3)
  52.                   Each of the elements in z must be positive.
  53.  
  54.      FINIT        (input) REAL
  55.                   The value of f at 0. It is more accurate than the one
  56.                   evaluated inside this routine (if someone wants to do so).
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSLLLLAAAAEEEEDDDD6666((((3333FFFF))))                                                          SSSSLLLLAAAAEEEEDDDD6666((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      TAU          (output) REAL
  75.                   The root of the equation f(x).
  76.  
  77.      INFO         (output) INTEGER
  78.                   = 0: successful exit
  79.                   > 0: if INFO = 1, failure to converge
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.